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Abstract. We investigate the evolution of pulsars in the (P, P) diagram. We first present analytical formulae to 
follow the evolution of a pulsar using simple exponential models for magnetic field decay and alignment. We then 
compare these evolutionary tracks with detailed model calculations using ohmic decay of crustal neutron star 
magnetic fields. We find that, after an initial phase with a small braking index, n, pulsars evolve with enhanced 
torque decay (n S> 3) for about 1 Myr. The long term evolution depends on the impurity parameter of the crust. 
If impurities are negligible in older isolated pulsars we expect their true age to be approximately equal to their 
observed characteristic age, r — P/{2P). It is not possible from data to constrain model parameters of the neutron 
star crust. 

Key words. Stars: neutron - Pulsars: general 



1. Introduction 

One of the key issues in pulsar physics is the question of 
torque decay. It is important to determine to what extent 
the braking torque decays in order to follow the dynamical 
evolution and calculate the true age of pulsars. The time 
dependent behaviour of the torque is closely linked to the 
important question of magnetic field decay and alignment 
of the magnetic and rotational axes of a neutron star. 

In Fig. 1 we have plotted the rotational period, P, ver- 
sus its derivative, P, for 947 observed pulsars (F. Camilo, 
private communications). The pulsars evolve to the right 
in this diagram as they loose rotational energy. To date, 
only four pulsars (see Lyne et al. 1996 and references 
therein) have a measured value of P (and the so-called 
braking index, n = UCt/Ct^, where 17 = 27r/P) which al- 
lows for a determination of the orientation of the evolu- 
tionary track in the (P, P) diagram for the given pulsar 
at its present age. However, these four pulsars are all very 
young (t = P/2P ~ 2 — 15 kyr) - otherwise P is too dif- 
ficult to measure. Thus we cannot learn much about long 
term torque decay from these observations. 
Here we will first derive analytical evolutionary tracks by 
considering simple B-field decay and alignment. Both of 
these effects will give rise to enhanced torque decay and 
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deviation from evolution along a straight line in the (P, P) 
diagram. Thereafter we shall present evolutionary tracks 
calculated from ohmic decay of crustal neutron star mag- 
netic fields. Such calculations of the neutron star crust 
were introduced by Sang & Chanmugam (1987) and Urpin 
& Muslimov (1992) who also included cooling models. 

The millisecond pulsars (P ^ 0.1 sec and P ^ lO"^'^) 
are seen to form a separate population in the (P, P) di- 
agram shown in Fig. 1 - presumably as a result of their 
binary origin (Alpar et al. 1982; van den Heuvel 1984). 
In this paper we shall only concentrate on the evolu- 
tion of isolated pulsars. For calculations of accretion- 
induced magnetic field decay of a binary neutron star 
(the millisecond pulsar progenitor) see e.g. Geppert & 
Urpin (1994), Urpin, Geppert & Konenkov (1997), Konar 
& Bhattacharya (1997; 1999). 

We outline the dynamics of a rotating neutron star in 
Sect. 2. In Sect. 3 we discuss a simple combined model with 
exponential magnetic field decay and alignment. In this 
section we also discuss the observed data and evidence for 
enhanced torque decay, as well as the role of the inclination 
angle with respect to the braking torque. A detailed model 
of the crustal physics of a cooling neutron star is used in 
Sect. 4 to calculate more realistic evolutionary tracks in 
the (P, P) diagram. Finally, the conclusions are briefly 
summarized in Sect. 5. 
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Fig. 1. Distribution of 947 observed pulsars in the {P,P) 
diagram. Binary pulsars are marked with a circle. The four 
pulsars with measured braking index evolve along the di- 
rection of the arrows. The solid line is an evolutionary 
track with n = 3. The dotted and dashed lines represent 
constant characteristic age and constant magnetic field, 
respectively. When radio pulsars cross the "death line" 
their emission is terminated as a result of the low elec- 
trostatic potential drop across their polar cap. In the top 
right corner is plotted observations of soft 7-ray repeaters 
(SGR) and anomalous X-ray pulsars (AXP). 

2. The rotating neutron star 

In the following section we shall concentrate on the "slow" 
(non-recycled) pulsars whose evolution is not "polluted" 
by interaction with a companion star. 
Consider a neutron star with moment of inertia, / and 
angular frequency, fl. Its rotational angular momentum 
and spin period is given by: 

L = m P = 2TT/fl (1) 
The rotational energy and its loss rate are given by: 

E,ot = ^/f^^ ^rot = (2) 

The braking torque acting on the neutron star is thus: 

N = ^ = m = E,,,/n (3) 

at 

(if we assume that / is constant in time). 

For an isolated neutron star E^q^ < 0, and hence N < 0. 



The change in magnitude of the torque as a function of 
time is given by: 



d\N\ 
dt 



-m 



(4) 



The term "torque decay" thus refers to a decrease in the 
magnitude of the torque. Hence, the needed condition is: 



dt 



< 



n > 



(5) 



where we have introduced the braking index, n defined as: 

(6) 



pp 



By introducing the characteristic age, r = P/{2P) we can 
get an idea of how | A''| decreases with the age of a pulsar: 

N = m = -T^I^ or |N| oc ^ (7) 



2.1. Torque decay and braking index in the 
magnetic-dipole model 

In the magnetic-dipole model (Pacini 1967, 1968; Ostriker 
& Gunn 1969) the spin-down energy is carried away by 
magnetic dipole radiation: 



3c3 



m = -- 



3c3 



(8) 



where R is the radius of the neutron star, B is the strength 
of the dipole component of the magnetic field at its equa- 
torial surface, m is the magnetic dipole moment and a 
is the inclination of the magnetic axis with respect to the 
rotation axis. Hence, in this model the braking torque (i.e. 
the radiation-reaction torque) transmitted to the neutron 
star by the magnetic field is proportional to B-^ sin^a. 
Prom the above equation one obtains: 



PP 



Stt^ 
3c37 



B^ sin^a 



(9) 



where we have assumed k i 



\c^I/ttR^ is a constant. 



3. Evolution with .B-field decay and alignment 

Soon after the discovery of pulsars Ostriker & Gunn (1969) 
suggested that the magnetic field of a neutron star should 
decay exponentially on a time-scale T£) ~ 4 Myr, because 
of ohmic dissipation of the supporting currents. The same 
authors also presented observational evidence in favour 
of field decay (Gunn & Ostriker 1970). However, it was 
soon argued by others that the ohmic decay cannot be 
important since the interior of a neutron star is likely to 
be superconducting (Baym, Pethick & Pines 1969). Ever 
since the early days of pulsar astronomy this important 
question of magnetic field decay has continued to be rather 
controversial, see e.g. Bhattacharya & Srinivasan (1995) 
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for a review on this matter. 

Let us for a moment assume an exponential field decaj^: 

B(t) = Bo e^*/^" (10) 

where tu is the time-scale of decay and Bq is the strength 
of the surface dipole magnetic field at time to — 0- 

Over the years, many authors have argued for 
or against observational evidence that the magnetic 
axis aligns with the rotation axis: Proszynski (1979), 
Candy & Blair (1983,1986), Lyne & Manchester (1988), 
McKinnon (1993), Gould (1994) and Gil & Han (1996). 
Beskin et al. (1988) even argue in favour of counter- 
alignment based on their magnetospheric theory. Recently 
Pandey & Prasad (1996) and Tauris & Manchester (1998) 
argue in favour of alignment. It is therefore important to 
investigate the evolution of the braking torque since this 
could be significantly influenced by the alignment process. 
The idea that the magnetic axis aligns with the rotation 
axis was first analyzed analytically by Jones (1976). He 
suggested a simple exponential decay of the inclination 
angle, a of the form: 



sinQ;(t) = sinofQ e 



(11) 



where ta is the alignment constant and ao is the initial 
inclination angle at the time to = 0- 

Inserting the two above expressions into Eq. (9) and 
integrating yields the period evolution of a pulsar in the 
general case including both exponential magnetic field de- 
cay and alignment: 

P'(0 =nofi5(l-e-2*/^")+ Po' (12) 

where Hq = Po-fo = {Bo sinao/fc)^. We have introduced 
an effective reduced time scale, td 
symmetry of Eqs (10) and (11). 
Using Eqs (3), (9), (12) and 17 = -27rF/p2 gives the time 
dependence of the braking torque: 



due to the 



N{t) 



-27r/Ho e 



-2t/fD 



[Hofz?(l-e-2Vr«) + p2]3/2 
A simple differentiation of Eq. (8) yields: 
^^^^ ^ _AR^Bit)B(t)suMt) ^3(,) 

oC 1 

AR^B'^it) sin a{t) cos a{t)a{t) 3 



3c3/ 

2R^B^{t) sin^ a{t) 



n%t) 



n^{t)n{t) 



and hence we have an expression for the braking index: 
3c^IB{t) 



n{t) 



R(^B^t) siii^ a{t) n^{t) 

3c^ I cot a{t) a (t) 
R^B'^jt) sin^ a{t) Q^jt) 



(15) 



^ This equation follows from Eq. (A. 8) in the appendix. 
From Eq. (A. 9) one finds td 4 Myr, as suggested by Ostriker 
& Gunn (1969), if L ~ 10** cm and a ~ 10^^ s'V However, this 
assumed value of the electrical conductivity, ( cr) is too small. 



or simply: 



n{t) = 3 - F{t) 



(16) 



(17) 



where 

Fit) = [3c'^I]/[R^B^{t)sm''a{t) n''{t)] 
and 

B{t) sin a{t) for B-decay and alignment 

s = < r ^.r--^' m 



B{t) for B-decay only 

sina(t) for alignment only 

Bo sinao for const. B and no alignment 



From the above equations we immediately obtain the al- 
lowed values of n{t), depending on the effect acting on it 
(cf. Table 9.1 of Manchester & Taylor 1977): 



> 3 for i? < and/or a < 
n{t) <( < 3 for (i > and B 

= 3 for B = Q and a = 



(19) 



where d > corresponds to the case of counter-alignment. 

3.1. Constant magnetic field and no alignment 

A Taylor expansion of Eqs (12), (13) and (15) for tjj 00 
or/and t^i ^ 00 yields the period, torque and braking 
index in case of a constant magnetic field, no alignment 
or both. For example, if we assume no field decay and no 
alignment (r^j 00 and ta 00) we obtain: 



P2(t) ^2 Hoi + Po' 



Nit) = 



-2ttII1, 







(2Hoi + P2)3/2 



(20) 
(21) 
(22) 



(13) 3.2. Results 

In the top panels of Fig. 2 we have plotted evolutionary 
tracks in the {P, P) diagram for various values of td and 
Ho- In the left panel we have plotted evolutionary tracks 
calculated without enhanced torque decay. The right panel 
shows evolutionary tracks including enhanced torque de- 
cay. For each initial starting point three curves are plotted 
corresponding to an increasing time scale, (from left 

(14) to right). We have assumed A; = 3.2 x 10^^ g2 cm 



s 2 . 



The small dots represent data of observed pulsars. It is 
seen that many evolutionary tracks fitting the data can 
be obtained from different choices of Ho and 775 . The thin 
lines in the top left panel are isochrones of 0, 1, 5, 10, 50, 
100, 500 kyr, 1, 5, 10, 50 and 80 Myr, respectively. 

To illustrate the difference between an evolution with- 
out {td — 00) and with {td = 10 Myr) enhanced torque 
decay (due to magnetic field decay, a multipole or align- 
ment) in Fig. 2 we have also plotted P, n, t and N, cal- 
culated from the equations derived earlier as a function of 
the true age, t. 
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Fig. 3. Dependence of the pulsar slow-down on the brak- 
ing index, n. The spin period, P is plotted as a function 
of true age, t for n = 2, n = 3 and n > 3 {td = 10 Myr). 
For n = 3 the slope of the curve is given by: P = P/{2t) 
and hence t = t (where _P » Pq = 10-0 ms). In Tabic 1 
is given the relation between true age, t and characteristic 
age, T for the different values of n. 

Table 1. Values of true age (t), characteristic age (t), 
spin period (P) and braking index n{t) for three different 
models of evolution (birth parameters: Pq = 10 ms and 
logPo = -12.5). 
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The evolution with enhanced torque decay (solid line, 
n > 3) is clearly seen to result in a strong decay of the 
braking torque for t > td and hence a strong decrease 
of P (P 0). Therefore a constant value of P is ap- 
proached as the radio emission process becomes inactive 
after 10-50 Myr (while the pulsar crosses the so-called 
death line). The strong decrease in P has an important 
influence on the braking index which accordingly increases 
sharply when t ^ fjj (see Eq. 6, P < 0). This is shown in 
the central right panel. 

If pulsars evolve without magnetic field decay or align- 
ment (dashed line, n = 3) they will attain much larger ages 
before they cross the death line. For these pulsars the true 
age, t will always be equal to the observed characteristic 



age, T = P/2P (P ^ Pq) - whereas if enhanced torque 
decay is important (n > 3) the true age of pulsars will 
be smaller than r (in this case by more than an order of 
magnitude for r > 20 Myr). 

In Fig. 3 we have demonstrated this effect graphically. 
The tangent to the n = 3 curve at a given age, t will al- 
ways intersect the x-axis (P = 0) at = —t. Thus the 
slope of this tangent is P{t) = P{t)/2t, or equivalently 
t = P(t) /2P{t) which is equal by definition to the char- 
acteristic age, T. The tangent to the n > 3 curve is seen 
to be less steep and therefore P{t) < P{t)/2t, or t < t. 
Finally, in the case of relaxed torque decay (0 < n < 3, 
e.g. radial deformation of B or a decreasing moment of 
inertia in young pulsars) we have t > t. 

For population statistics and evolution of pulsars it 

is therefore essential to know whether or not enhanced 
torque decay occurs in order to determine the true ages of 
the observed pulsars. 

It should be noted, that strictly speaking the braking in- 
dex, n is defined from the power-law: f2 oc — f2". However, 
in general n is a function of time: O = A{t) and since 
there is no analytical solution to this power-law equation, 
we have retained the definition of n as given in Eq. (6). 



3.3. Observational evidence on torque decay 7 

Many attempts have been made to find observational ev- 
idence for magnetic field decay. The key problem (as out- 
lined above) is that we do not know the true age of the pul- 
sars observed. Furthermore, many parameters connected 
to pulsars (e.g. r and B) are given as functions of the two 
fundamental observed parameters: P and P. There are 
many different ways to illustrate the evolution of pulsars. 
Here, we introduce another variant of the (P, P) diagram 
and plot the rotational energy of a pulsar as a function of 
its characteristic age. This is shown in Fig. 4. The dashed 
lines represent evolution with n = 3 (constant B-field and 
magnetic inclination angle) and the solid lines are cal- 
culated with enhanced torque decay (n > 3, fi) = 10 
Myr). In each case we considered two different sets of 
initial birth parameters for the pulsar: (Pq = 10.0 ms, 
logP = -13.5 ^ Bo = 5.69 x lO^^ G) and (Pq = 31.6 ms, 
logP = -12.0 Bo = 5.69 X lO^^ q) jf pulsars evolve 
without enhanced torque decay they should move along 
the dashed lines. If enhanced torque decay is important 
their evolution should bend off along the solid lines. 

We have shaded two regions in the diagram (A and B) 
where there is a deficit of pulsars. If pulsars evolve with 
n = 3, the region A should feed pulsars into the relatively 
dense populated area near the arrowhead. Similarly, re- 
gion B should be the endpoint of pulsars coming from the 
populated area along its arrow base. So why are there so 
few pulsars in region A and B ? A simple reason for the 
few pulsars in region A could be due to the logarithmic 
scale for the characteristic age. Hence, the pulsars accu- 
mulate with time near the arrowhead before crossing the 
death line. However, if this is the correct explanation it 
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Fig. 4. The rotational energy of observed pulsars plotted as a function of their characteristic age. The solid and dashed 
lines are evolutionary tracks with and without enhanced torque decay, respectively. We assumed a pulsar moment of 
inertia, / = 10^^ g cm^. The millisecond pulsar population is seen in the upper right corner. See text for discussion. 



seems strange why there should not be a similar effect 
below this region, e.g. between and parallel to the dashed 
lines. Furthermore, the deficit of pulsars in region B is even 
more strange in this respect (since pulsars are then sup- 
posed to be accumulated in this region), unless the naively 
drawn death line in reality bends into this region. In both 
cases we notice that a natural explanation, for the deficit 
of pulsars in the shaded regions, is the effect of enhanced 
torque decay and subsequent bending of the evolutionary 
tracks. We will let the reader draw his own conclusion. 

3.3.1. Torque decay and the magnetic inclination angle 

In the Goldreich- Julian (1969) model of a non-vacuum 
magnetosphere of a pulsar, the slow-down of rotation is 
caused by outflow of plasma: 

^™t = -(^-^) ^ttRIc (23) 

where i?L = c/^ is the radius of the light cylinder and 
is the magnetic field strength at this location - for 



a dipole field: Bl oc or = B{R/Ri^f. This ex- 
pression for E^ot is of the same magnitude as in Eq. (8), 
but here the braking torque (N = E-cot/^) is indepen- 
dent of a. The dynamics of a rotating neutron star in 
this case is then similar to that of the magnetic-dipole 
model for a constant inclination angle. In Fig. 5 we have 
plotted the braking torque, N oc P/P^ vs. the observed 
inclination angle, Ofobs for a large number of pulsars. The 
derived, observed magnetic inclination angles were taken 
from the work by Tauris & Manchester (1998) and refer- 
ences therein. Though we cannot rule out the possibility 
that N is independent of a (a x^-test yields a ~ 15 % 
probability for constancy) we notice that N seems to be 
a slightly increasing function of a (although the scatter is 
large )^. However, we must also bear in mind that a itself 
is only a marginally decreasing function of characteristic 



This is in contrast to Bhattacharya (1989) who argued that 
the spin-down torque is independent of a. However, he consid- 
ered explicitly the dipole model and also assumed a constant 
magnetic field strength. This is why he plotted PP, rather than 
A'^ = itl cx P /P^ , as a function of a. 
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Fig. 5. A plot of P/P^ vs. the observed inclination angle, 
aobs- The top panel contains 316 pulsars obtained from 
Gould (1994) and the bottom panel shows 149 pulsars 
from Rankin (1993) - see Tauris & Manchester (1998) for 
details. The braking torque, N oc P / P"^ is seen to be a 
slightly increasing function of aobs- Standard error bars of 
the binned mean values are plotted. 

age (see Tauris & Manchester 1998). Hence, based on ob- 
servations the effect of a on iV is probably not a major 
concern. 

We have now demonstrated how simple analytic mod- 
els can account for the distribution of observed pulsars in 
the (P, P) diagram, although the behaviour of the evolu- 
tionary tracks is still an unsettled question. The models 
we have used so far are based on crude approximations 
(e.g. the exponential decay of the S-field). 
In the following section we will apply more realistic physics 
to our calculations and see how this will affect the calcu- 
lated evolutionary tracks. 

4. Ohmic decay of crustal neutron star B-fields 

The physics and relevant equations of crustal i?-field 
decay is given in the appendix. We refer to Konar & 
Bhattacharya (1997; 1999) and Konar (1997) and refer- 
ences therein for further details on the computations. Here 
we shall briefly summarize the assumptions necessary for 
the calculations we have performed in this papcr^. 

^ For an alternative scenario of neutron star magnetic field 
evolution caused by core field changes and crust movements, 
see e.g. Ruderman, Zhu & Chen (1998). 



We assume that the magnetic field has been gener- 
ated in the outer crust by some unspecified mechanism, 
e.g. by thermomagnetic effects (Blandford, Applegate & 
Hernquist 1983), during or shortly after the neutron star 
is formed. For the distribution of the initial underlying 
toroidal current we have chosen to use the configuration 
by Konar & Bhattacharya (1997). 

The equation of state by Wiringa, Fiks & Fabrocini (1988) 
is matched to Negele & Vautherin (1973) and Baym, 
Pcthick & Sutherland (1971) in order to construct a neu- 
tron star with a mass of 1.3 M© (R — 11.1 km). 
The electrical conductivity in the neutron star crust is 
mainly a function of mass density and temperature. The 
conductive properties of the solid crust are determined by 
scattering of electrons on phonons and lattice impurities 
(Yakovlev & Urpin 1980 and Itoh et al. 1993). In the outer 
thin liquid layer near the surface it is determined by scat- 
tering of electrons on ions (e.g. Urpin & Muslimov 1992). 
We used the models by van Riper (1991a, 1991b) for 
standard cooling of a neutron star (modified URCA pro- 
cess in the core and neutrino pair bremsstrahlung in the 
crust). We adopted the profile of Gudmundsson, Pethick 
& Epstein 1982) to estimate the temperature in the 
outer/surface regions {p < 10^" g cm^"^) above the ap- 
proximately isothermal inner crust/core region. 

We have performed calculations of the ohmic decay of 
the crustal magnetic field of an isolated pulsar for an age 
up to 100 Myr. In Fig. 6 we have plotted the evolution in 
the (P, P) diagram. Here we assumed that the pulsar is 
born with an initial period, Pq — 10 ms and initial period 
derivative, Pq = 10~^^, corresponding to an initial surface 
B-field strength of Bo = 10^^ Gauss. The three different 
main curves are the results of calculations using different 
initial locations (densities) of the current configuration in 
the crust. The curves from top to bottom correspond to 
logpo = 13, 12 and 11 (g cm~^), respectively. The branch- 
ing of the curves near their endpoints of evolution is caused 
by different assumed values of the impurity parameter: 

Q = n^iZ ~ Z^f (24) 

"total ^ 
I 

where ntotai is the total ion density, ri; is the density of im- 
purity species i with charge Zi, and Z is the ionic charge in 
the pure background lattice (Yakovlev & Urpin 1980). The 
four branches correspond to: Q = 0, 0.01, 0.05 and 0.10 
(top to bottom, respectively). As the neutron stars cools 
sufficiently, increasing values of Q will decrease the electri- 
cal conductivity of the crust. In Fig. 7 we demonstrate how 
a single model with fixed parameters (po = 10^^ g cm"'^ 
and Q = 0) can cover the entire population of observed 
non-recycled pulsars by choosing different initial values of 
Pq (here corresponding to logSo = 13.5, 13.0 and 12.5, 
keeping Pq — 10 ms fixed). 

The evolutionary tracks obtained from our calculations 
of neutron star crustal field decay are also seen to cover 
the observed distribution of pulsars in the (P, P) diagram 
quite well. However, from Figs. (6) and (7) it is also clear 
that the location of each track depends on the underlying 
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Fig. 6. Evolutionary tracks in the (P, P) diagram assum- 
ing Po = 10 ms and Pq = 10"" (^o = 10" G) calculated 
from Po = 10" - 1013 g cm'^ and Q = - 0.10 ~ see 
text. Each track was followed for 100 Myr. The small dots 
represent data. 



parameters of the applied model. It is impossible to con- 
strain e.g. Po and Q from the observed data set. Whereas 
a spread in the initial magnetic field strengths. Bq can be 
expected, other parameters, like the impurity parameter, 
Q, should be more or less equal for all pulsars. 

The calculated evolutionary tracks, of pulsars with 
ohmic crustal magnetic field decay, are seen to initially 
evolve along straight (n = 3) lines in the (P, P) diagram 
(Figs. 6 and 7), then bending down (n > 3) for a while and 
returning again to an evolution along a straight (n = 3) 
line, if Q = 0. This is also shown in Fig. (8) where we have 
plotted the braking index, n as a function of true age. t. 
The evolution of n as a function of log t can roughly be de- 
scribed to be Gaussian-like for the first 10 Myr. Initially, 
n = 3 (by assumption), then follows a short interval with 
3 < n < 30 before it settles down again and evolves with 
n » 3. The more exact evolution in this temporary inter- 
val (0.01 ^ t ^ 1 Myr) depends on the initial location 
(depth) of the supporting currents, po and the details of 
the cooling process (exotic vs. standard cooling). The rea- 
son for the flattening of the evolutionary track (returning 
to n = 3) at ages t > 1 Myr, is due to a saturation of 
the decay of the B-field. The current distribution, which 
is responsible for the B-field, migrates inward as a result 
of diffusion and enters the highly conducting parts of the 
neutron star. In that region rohmic > THubbie and hence 
the magnetic field will essentially be stable forever - it 
freezes out at a residual value (Konar 1997). 




0.01 0.1 



P [s] 

Fig. 7. Evolutionary tracks calculated for po = 10^^ g 
cm"'^ and Q = at different initial B-fields for Pq = 10 
ms. 
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Fig. 8. The braking index as a function of true age calcu- 
lated for a pulsar with crustal magnetic field decay. We as- 
sumed Bo = 10"G (Po = 10 ms), Q = 0.0 and log/3o=ll, 
12 and 13 (g cm"^). The star indicates an average value 
obtained from Johnston & Galloway (1999) - see text for 
discussion. 



Since the x-axis in Fig. 8 is logarithmic, a pulsar with 
Q — spends the majority of its life as an active pulsar 
evolving with n « 3. Hence, near the end of the evolution- 
ary tracks {t « 100 Myr) the characteristic age, t is only 
~ 3% larger than the true age, t. 

The star in the centre of the figure is the mean value of n 
for the eight pulsars with the smallest error bars in Table 1 
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Fig. 9. The braking index as a function of true age cal- 
culated for a pulsar with crustal magnetic field decay. We 
assumed Bq = 10"G (Pq = 10 ms), po = lO^^ g ^^^-3 
and (5 = 0, 0.01, 0.05 and 0.10. See text for discussion. 



enhanced torque decay (n 3> 3). This can explain the dis- 
tribution of observed pulsars in the {P, P) diagram quite 
well. The subsequent evolution depends on the impurities 
in the crustal crystalline lattice of old isolated pulsars. If 
Q — the older pulsars (t > 10 Myr) returns to long term 
n ~ 3 evolution and future observations of the temper- 
ature of nearby, relatively old pulsars (e.g. PSR J0108- 
1431, Tauris et al. 1994) may be able to verify or reject 
our hypothesis that t « r for older pulsars. 

Further investigations of the present pulsar data in the 
(P, P) diagram and statistical population analyses are en- 
couraged in order to help verify our finding of enhanced 
torque decay at some stage during the pulsar evolution. 
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of Johnston & Galloway (1999). Note, that they quote 
characteristic pulsar ages, t and we have not corrected for 
the unknown reduction in the true ages, t. Although there 
are large uncertainties involved, it is interesting to notice 
that our calculations predict values of the braking index, 
n within the range of their estimates. 

It has recently been suggested that the impurity pa- 
rameter is probably quite large in a heated crust of an 
accreting binary neutron star as a result of nuclear pro- 
cesses (Brown & Bildsten 1998). However, it is uncertain 
whether or not Q is important for isolated pulsars. If there 
are significant impurities {Q > 0) in the crustal crystalline 
lattice of isolated pulsars, the braking index will increase 
again after 10 Myr. This is demonstrated in Fig. 9 (see 
also the final branching of the curves in Fig. 6). Hence, 
a determination of the braking index for such old pulsars 
could in principle yield the value of Q. Unfortunatly such 
direct measurements of n for old pulsars are not possible 
(see e.g. Johnston & Galloway 1999). 

5. Conclusion 

The effects of long term torq\.ic decay are presently not 
known very well. The only pulsars with a measurable brak- 
ing index are quite young. If more young pulsars turn out 
to have 2 ^ n < 3 it seems clear that n must be increasing 
with age. Otherwise, their evolution will not be consistent 
with the bulk of observed pulsars with ages 1 — 10 Myr. 
As we have argued in this paper, enhanced torque decay 
{n > 3) seems to play an important role for pulsars. The 
exponential models for B-field decay and alignment are 
obviously naive. However, they give a better description 
of the observed data than those models with n = 3. 

Using a detailed model of the crustal physics of a neu- 
tron star, we have demonstrated that pulsars may evolve 
for a limited interval of time (from 0.01 — 1 Myr) with 



Appendix A: Physics of the S-field decay 

The induction equation is given by (Jackson, 1975): 

V = -— Vx(-x VxB) + Vx(FxB) (A.l) 

at iTT a 

where V is the velocity of material movement and a is 
the electrical conductivity of the medium. We choose the 

vector potential to be of the form A = (0, 0, A^) to ensure 
a poloidal geometry for B = V x A. In. particular: 



g{r, t) sm.9 



(A.2) 



where g{r,t) is the Stokes' stream function. In the lowest 
order of multipole, the dipolar form of the B-field is: 



B{r,e) = V X ( 



g{r, t) sva.6 



2cos0g{r,t) ^ sm9dg{r,t)^ 
The magnitude of the field is then given by: 



B{r,e) 



A cos^ eg^{r,t) sin^ dg{r,t) 



+ 



-(- 



1/2 



(A.3) 



(A.4) 



r"* dr 
Therefore, at the pole {9 = Q, r = R) the magnitude is: 

2g{R,t) 



B{R,0) = 



P2 



(A.5) 



which is simply proportional to the value of the Stokes' 
function there. The underlying current distribution corre- 
sponding to the above field is given by: 



J = —V X B 
An 

c sin 6* 
47r r 



d^gir,t) 2g{r,t) 



Qj.2 



(A.6) 
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A.l. Pure ohmic diffusion 
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In order to understand the ohmi c dis sipation of the field 
strength, let us consider equation [ A.l | without the second 
term on the right hand side. The second term corresponds 
to the convcctive transport, which is relevant for material 
movement through the crust of an accreting neutron star. 
Without this term equation [|A.l| takes the following form: 

no „2 1 

_.-_Vx(-xVxB) (A,7) 

For the uniform conductivity case, the above equation 
takes the form of a pure diffusion equation (by virtue of 

the divergence-free condition for the magnetic fields): 

dB _ 

~dt ^ ^ Ana 

The diffusion constant for the above equation is 
and the time-scale characteristic of the process is: 
47rcrL2 



Tdiff 



(A.8) 

V47ra 



(A.9) 



where L is the length-scale associated with the underlying 
current distribution supporting the field. 



A. 2. The field evolution equation 



We use the form of B{r,6,4>) given by equation |A.3| to 
cast equation | A.7] in terms of the Stokes' stream function. 
I. The left hand side: 
dB _ dVxA 
'dt ~ di 



V X 



dt 



g(r,0) sint 



(A.IO) 



II. The right hand side: 



V X 



V X 



V X 



-V X B 

a 



1^ 
-V X 

a 



2 cos 6* g{r, t) 



1 sine 
cr r 



dgjr.t) 
dr 



d^g{r,t) 2g{r,t) 



(A.ll) 



Incor pora tion of the expressions [ A.IC | and [ A.ll | in equa- 
tion IA.TI then leads to: 



dg{r,t) 
dt 



c 



d^gjr.t) 2g{r,t) 



(A.12) 



The results presented in this paper are based on numer- 



ical solutions of the equation [A.12|. To solve this equa- 
tion we used a modified Crank-Nicholson scheme of dif- 
ferencing (Press et al. 1992) and solved the relevant tri- 
diagonal matrices. We used the standard boundary con- 
ditions, e.g. introduced by Geppert & Urpin (1994). We 
refer to Konar (1997) for an extended discussion on the 
details of obtaining a stable mathematical solution, and 
the physics of the boundary conditions. 
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